THE LAURENT PHENOMENON 



SERGEY FOMIN AND ANDREI ZELEVINSKY 



Abstract. A composition of birational maps given by Laurent polynomi- 
als need not be given by Laurent polynomials; however, sometimes — quite 
unexpectedly — it does. We suggest a unified treatment of this phenomenon, 
which covers a large class of applications. In particular, we settle in the affir- 
mative a conjecture of D. Gale and R. Robinson on integrality of generalized 
Somos sequences, and prove the Laurent property for several multidimensional 
recurrences, confirming conjectures by J. Propp, N. Elkies, and M. Klcbcr. 
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1. Introduction 

In this paper, we suggest a unified explanation for a number of instances in which 
certain recursively defined rational functions prove, unexpectedly, to be Laurent 
polynomials. We begin by presenting several instances of this Laurent phenomenon 
established in the paper. 

Example 1.1. (The cube recurrence) Consider a 3-dimensional array 

(Vijk ■ {i,j,k) G H) 
whose elements satisfy the recurrence 

ctyi-i,j,kVi,i-i,k-\ + 0Vi,j-i,kyi-i,j,k-i + jy%,j,k-iyi-i,j-i,k 



(l.i) 



Vi,j,k 



yi-l,j-l,k-l 

Here TL can be any non-empty subset of Z 3 satisfying the following conditions: 

(1.2) if (i,j,k) G H, then (i',j',k') G H whenever i < < j',k< k'\ 

(1.3) for any (i',j',k') e H, the set {(i,j,k) eH:i<i',j< j',k < k'} is finite. 
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Theorem 1.2. Let H iui t = {(a,b,c) G U : (a - 1, b - 1, c - 1) ^ H}. For every 
(i,j,k) G 7i, i/ie entry y%^,k is a Laurent polynomial with coefficients in "L[a,(3,^\ 
in the initial entries y a .b,c, for (a, 6, c) G -f/init- 

The cube recurrence (with a = /3 = 7 = l) was introduced by James Propp JlOj] , 
who was also the one to conjecture Laurentness in the case when Ti C Z 3 is given 
by the condition i + j + fc > 0; in this case -ffmit consists of all (a, b,c) £ Ti such 
that a + b + c £ {0,1, 2}. Another natural choice of Ti was suggested by Michael 
Kleber: Ti = Z> , in which case i?i n it = {{a, b, c) G Z> : a6c = 0}. 

Example 1.3. (The Gale-Robinson sequence) Let p, q, and r be distinct positive 
integers, let n = p + q + r, and let the sequence yo, yx, • • • satisfy the recurrence 

,x ay k+p y k + n ~p + (3y k + q yk+n- q + lVk+ r yk+n-r 
(1.4) y fc+n = ^ ^ 



yk 

David Gale and Raphael Robinson conjectured (see and || E15]) that every 
term of such a sequence is an integer provided yo = ■■ ■ = y n -i = 1 and a, (3, 7 are 



positive integers. Using Theorem 1.2, we prove the following stronger statement. 



Theorem 1.4. As a function of the initial terms yo, . . . , y n -\, every term of the 
Gale-Robinson sequence is a Laurent polynomial with coefficients in Z[a,f3,j]. 

We note that the special cas e a — f3 = j = 1, p = 1, q = 2, r = 3, n = 6 (resp., 
r = 4, n — 7) of the recurrence ( |l.4|) is the Somos-6 (resp., Somos-7) recurrence JtJ. 

Example 1.5. (Octahedron recurrence) Consider the 3-dimensional recurrence 

r\ _ a yi+i,j,k-iyi-i,j,k-i + (3yi,j+i,k-iyij-i,k-i 

I.J--OJ Ui,j,k — 

for an array {yijk)u,j,k)EH whose indexing set Ti is contained in the lattice 
(1.6) L = {(i,j,k) G Z 3 : i + j + k = mod 2} 



and satisfies the following analogues of conditions fll.2| )-(1.3): 

(1.7) if (i,j,k) £ Ti, then {i',j',k') G Ti whenever \i' -i\ + \f - j\ < k' - k; 

(1.8) for any (i',j',k') G Ti, the set {(i,j,k) G Ti : \i' - i\ + \f - j\ < k' - k} 
is finite. 



Theorem 1.6. Let H init = {(a,b,c) £ Ti : (a,b,c— 2) ^ Ti}. For every (i,j,k) G 
Ti, the entry yij lk is a Laurent polynomial with coefficients in 2\a,(3] in the initial 
entries y a .b,c, for (a,b,c) G -ff in it- 

The octahedron recurrence on the half-lattice 

(1.9) Ti = {(i,j,k) £L : fc>0} 

was studied by W. H. Mills, D. P. Robbins, and H. Rumsey in their pioneering 
work H on the Alternating Sign Matrix Conjecture (cf. Q and |l^, Section 10] for 



further references); in particular, they proved the special case of Theorem 1.6 for 
this choice of Ti. 
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Example 1.7. (Two-term version of the Gale-Robinson sequence) Let p, q, and n 
be positive integers such that p < q < n/2, and let the sequence yo,yi, ■ ■ . satisfy 
the recurrence 

(1-10) Vk+n = " 1 -• 

Vk 



Using Theorem 1.6, one can prove that this sequence also exhibits the Laurent 
phenomenon. 

Theorem 1.8. As a function of the initial terms yo, . . . , y n -i> every term y m is 
a Laurent polynomial with coefficients in Z[a, 0\. 

We note that in the special case a = (3 = 1, p = 1, q = 2, n = 5 (resp., n = 4), 
( |1.1C| ) becomes the Somos-5 (resp., Somos-4) recurrence 0. 

The last example of the Laurent phenomenon presented in this section is of a 
somewhat different kind; it is inspired by 

Example 1.9. Let n > 3 be an integer, and consider a quadratic form 



P(xi, . . . ,X n ) = xj-\ \-X% +^ CtijXiQ 



i<j 

Define the rational transformations F\ , . . . , F n by 

P\ 

(1.11) Fi : (xi, . .., x n ) ' ^ (xi, . . . , Xi-i, — ^^,Xi + i, . . . , x„). 

Theorem 1.10. For any sequence of indices i±, . . . , i m , the composition map G = 
Fi t o • ■ ■ o Fi m is given by 

G : x = (x\, ...,x n )^> (Gi(x), . . . , G n (x)), 

where G\, . . . , G n are Laurent polynomials with coefficients in Z[(Xij : i < j]. 

This paper is an outgrowth of where we initiated the study of a new class 
of commutative algebras, called cluster algebras, and established the Laurent phe- 
nomenon in that context. Here we prove the theorems stated above, along with a 
number of related results, using an approach inspired by The first step is to 
reformulate the problem in terms of generalized exchange patterns (cf. Q Defini- 
tion 2.1]), which consist of clusters and exchanges among them. The clusters are 
distinguished finite sets of variables, each of the same cardinality n. An exchange 
operation on a cluster x replaces a variable x 6 x by a new variable x' = — , where 
P is a polynomial in the n — 1 variables x — {x}. Each of the above theorems can be 
restated as saying that any member of the cluster obtained from an initial cluster 
Xq by a particular sequence of exchanges is a Laurent polynomial in the variables 



from Xq. Theorem 1.1C is explicitly stated in this way; in the rest of examples 



above, the rephrasing is less straightforward. 

Our main technical tool is "The Caterpillar Lemma" (Theorem 2.1), which es- 
tablishes the Laurent phenomenon for a particular class of exchange patterns (see 
Figure [l]). This is a modification of the namesake statement || Theorem 3.2], 
and its proof closely follows the argument in (We note that none of the two 
statements is a formal consequence of another.) 
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In most applications, including Theorems 1.2 and 1.6 above, the "caterpillar" 
patterns to which Theorem 2.1 applies, are not manifestly present within the origi- 
nal setup. Thus, we first complete it by creating additional clusters and exchanges, 
and then apply the Caterpillar Lemma. 

The paper is organized as follows. The Caterpillar Lemma is proved in Sec- 
tion 2| Subsequent sections contain its applications. In particular, Theorems 1.2 



1.4 



tion 



1.6, and 1.8 are proved in Section ^, while Theorem 1.10 is proved in Sec 
5[~Othcr instances of the Laurent phenomenon treated in this paper include 
generalizations of each of the following: Somos-4 sequences (Example 3.3), Elkies's 
"knight recurrence" (Example 4.1), frieze patterns (Example |4.3| ) and number walls 
(Example O) . 

We conjecture that in all instances of the Laurent phenomenon established in this 
paper, the Laurent polynomials in question have nonnegative integer coefficients. 
In other contexts, similar nonnegativity conjectures were made earlier in @, |^, 

Acknowledgments. We thank Jim Propp for introducing us to a number of 
beautiful examples of the Laurent phenomenon, and for very helpful comments on 
the first draft of the paper. In particular, it was he who showed us how to deduce 
Theorem 1.8 from Theorem 1.6. 

This paper was completed during our stay at the Isaac Newton Institute for 
Mathematical Sciences (Cambridge, UK), whose support and hospitality are grate- 
fully acknowledged. 



2. The Caterpillar Lemma 

Let us fix an integer n > 2, and let T be a tree whose edges are labeled by the 
elements of the set [n] = {1, 2, . . . , n}, so that the edges emanating from each vertex 
receive different labels. By a common abuse of notation, we will sometimes denote 

by T the set of the graph's vertices. We will write t — t' if vertices t, t' S T are 

joined by an edge labeled by k. 

From now on, let A be a unique factorization domain (the ring of integers Z 
or a suitable polynomial ring would suffice for most applications). Assume that 
a nonzero polynomial P £ A[x±, . . . , x n ], not depending on x k , is associated with 

every edge t — ^— t' in T. We will write t — — t' or t — — t' , and call P the 

exchange polynomial associated with the given edge. The entire collection of these 
polynomials is called a generalized exchange pattern on T. (In 0, we introduced a 
much narrower notion of an exchange pattern; hence the terminology.) 

We fix a root vertex to € T, and introduce the initial cluster x(to) of n in- 
dependent variables xi(io), ■ ■ ■ , x n {fo). To each vertex t 6 T, we then associate a 
cluster x(i) consisting of n elements xi(t), . . . , x n (t) of the field of rational functions 
A(xi(to) 1 • ■ ■ , x n (to)). The elements Xi(t) are uniquely determined by the following 

exchange relations, for every edge t — — t': 

(2.1) Xi(t) = Xi(t') for any i ^ k; 

(2.2) x k (t)x k (t r )=P( X {t)). 

(One can recursively compute the ^(i)'s, moving away from the root. Since the 
exchange polynomial P does not depend on Xk, the exchange relation ( |2.2| ) does 
not change if we apply it in the opposite direction.) 
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We next introduce a special class of "caterpillar" patterns, and state conditions 
on their exchange polynomials that will imply Laurentness. 

For to > 1, let T n . m be the tree of the form shown in Figure B. 



to ihcad 

7T7V7V7V7V7Y 



Figure 1. The "caterpillar" tree T njTO , for n — 4, to = 8 

The tree T„ jm has m vertices of degree n in its "spine" and m(n — 2) + 2 vertices 
of degree 1. We label every edge of the tree by an element of [n], so that the n 
edges emanating from each vertex on the spine receive different labels. We let the 
root to be a vertex in T nm that does not belong to the spine but is connected to 
one of its ends. This gives rise to the orientation of the spine, with all the arrows 
pointing away from to (see Figure |l|). We assign a nonzero exchange polynomial 

P G A[xi, . . . , x n ] to every edge t t' of T„. m , thus obtaining an exchange 

pattern. 

For a rational function F — F{x, y, . . . ), we will denote by F\ x <_ g ( x ,y,...) the 
result of substituting g{x, y, . . . ) for x into F. To illustrate, if F{x, y) = xy, then 
FL^v = £. 

Theorem 2.1. (The Caterpillar Lemma) Assume that a generalized exchange pat- 
tern on T n ,m satisfies the following conditions: 

(2.3) For any edge • — — •, the polynomial P does not depend on Xk, and is not 
divisible by any Xi, i S [n]. 

(2.4) For any two edges • • — > •, the polynomials P and Qo — Q\ Xi =o 
are coprime elements of A[xi, . . . , x n ]. 

(2.5) For any three edges • -~- • • -~- • labeled we have 

L ■ Q ■ P = R\ , 

where b is a nonnegative integer, Qq — Q\ Xi =a > an d L * s a Laurent 
monomial whose coefficient lies in A and is coprime with P. 

Then each element Xi(t), for i £ [n], t € T nm , is a Laurent polynomial in 
x\(to), . . . , x n (to), with coefficients in A. 



(Note the orientation of edges in (2.4 )— (|2.5|) .) 



Proof. Our argument is essentially the same as in |g, Theorem 3.2]. For t € T„. m , 
let 



C(t)=A[x 1 (t) ±1 ,...,x n (t) ± ^ 
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denote the Laurent polynomial ring in the cluster x(t) with coefficients in A. We 
view each C{t) as a subring of the ambient field of rational functions A(x(t )). 

In this notation, our goal is to show that every cluster x(t) is contained in £(to). 
We abbreviate Co — C(to). Note that £ is a unique factorization domain, so 
any two elements x,y S Co have a well-defined greatest common divisor gcd(x,y) 
which is an element of Co defined up to a multiple from the group Cq of invertible 
elements in Co; the group Cq consists of Laurent monomials in x±(to), ■ ■ ■ ,x n (to) 
whose coefficient belongs to A x , the group of invertible elements of A. 

To prove that all x(i) are contained in Co, we proceed by induction on m, the 
size of the spine. The claim is trivial for m — 1, so let us assume that m > 2, and 
furthermore assume that our statement is true for all "caterpillars" with smaller 
spine. It is thus enough to prove that x(thcad) C Co , where thead is one of the 
vertices most distant from to (see Figure |l|) . 

We assume that the path from to to thcad starts with the following two edges: 

to -~- ti — ^— > t%. Let £3 S T„ m be the vertex such that t 2 — ^— t 3 . The following 

P Q R 

lemma plays a crucial role in our proof. 

Lemma 2.2. The clusters x(ii), x(t 2 ), and ^c(t 3 ) are contained in Co- Further- 
more, gcd{x l (t 3 ),x i (t 1 )) = gcd(Xj(t 2 ),Xi(ti)) = 1. 



Proof. The only element in the clusters x(ti), x(t2), and ^(t^) whose inclusion 
in Co is not immediate from ( |2.l| )-( |2~2] ) is Xi(t 3 ). To simplify the notation, let us 
denote x = x l (t ), y = Xj(t ) = Xj(h), z = Xi(t{) = x l (t 2 ), u = Xj(t 2 ) = Xj(t 3 ), 
and v = Xi(ts), so that these variables appear in the clusters at to, . . . , t^, as shown 
below: 



y,x z.y u,z v,u 

• • > • • 

P Q R 

to tl t2 t3 

Note that the variables Xk, for k ^ {i,j}, do not change as we move among the 
four clusters under consideration. The lemma is then restated as saying that 

(2.6) v e C ; 

(2.7) gcd(« J «) = l; 

(2.8) gcd(z,«) = l. 

Another notational convention will be based on the fact that each of the polynomials 
P, Q, R has a distinguished variable on which it depends, namely Xj for P and R, 
and Xi for Q. (In view of ( |2.3| ), P and R do not depend on Xi, while Q does 
not depend on Xj.) With this in mind, we will routinely write P, Q, and R as 
polynomials in one (distinguished) variable. For example, we rewrite the formula 



in (2.5) as 



(2.9) R [ ) = L(y)Q(0) b P(y), 



y 



where we denote L(y) = L\ Xj <— y . In the same spirit, the notation Q' , R' , etc., will 
refer to the partial derivatives with respect to the distinguished variable. 
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We will prove the statements (2.6), (2.7), and (2.8) one by one, in this order. 
We have: 



Q{z) 



R(u 



y 

'Q(z) 



R 



(T) 



R 



Since 



and 



R 



e Co 



R ( 9 V 1 ) L(y)Q(0) b P(y) 



(2.6) follows. 



We next prove (2.7). We have 



Q(z) _ Q(0) 



L(y)Q(0) b x e C , 



mod 



y y 

Since x and y are invertible in Co, we conclude that gcd(z, u) = gcd(P(y), Q(0)) = 1 
(using (H|)). 

It remains to prove (|2.8|). Let 



f{z) = R (Q& 



Then 



Working modz, we obtain: 



m-w) 



L(y)Q(0) b x. 



/(0) 



Hence 



v J 



no) 



Q'(Q) 



Q'(Q) 
v 



L{y)Q(0) b x mod z. 



Note that the right-hand side is a polynomial of degree 1 in a; whose coeffici ents 
are Laurent polynomials in the rest of the variables of the cluster x (tp). Thus ( [2.8|) 
follows from gcd (L(y)Q(0) b , P(y)) — 1, which is a consequence of (2.4)- (|2.5|) . □ 

We can now complete the proof of Theorem 2.1. We need to show that any 
variable X = Xfc(thead) belongs to £o- Since both t\ and £3 are closer to ihcad 
than to, we can use the inductive assumption to conclude that X belongs to both 
C(tx) and £(£3). Since X 6 £(fi), it follows from (2.1) that X can be written as 
X = f /xi(ti) a for some / G Co and a 6 Z> . On the other hand, since X e C(t 3 ), 
it follows from (2.1) and from the inclusion Xi{tz) 6 £0 provided by Lemma 
that X has the form X = g / Xj{t2) b Xi{t^) c for some g £ Co and some b, c G Z>o - 
The inclusion X e £0 now follows from the fact that, by the last statement in 
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Lemma 2.2, the denominators in the two obtained expressions for X are coprime 



in Cq. □ 
3. One-dimensional RECURRENCES 



In this section, we apply Theorem 24 to study the Laurent phenomenon for 
sequences yo,yx,--- given by recursions of the form 

(3.1) Vm+nVm = F(y m+ x, y m +n-l), 

where F € A[xi, . . . , x n -x]- 

For an integer m, let (m) denote the unique element of [n] = {1, . . . , n} satisfying 
m = (m) mod n. We define the polynomials Fx, . . . ,F n £ A[xi, . . . , x„] by 

(3.2) F m = F{x 

(m+l> , x {m+2) , ■ ■ ■ , X 

(m-l>); 

thus F m does not depend on the variable x m . We introduce the infinite "cyclic 
exchange pattern" 

(3-3) t t x t 2 -®— t 3 -J®— U , 

-f*(0> *"<1) -T{2> -T(3) 

and let the cluster at each point t m consist of the variables y m , . . . , y m + n -i , labeled 
within the cluster according to the rule y s = X/ s \ (t m ). Then equations ( |3.l| ) become 
the exchange relations associated with this pattern. 

To illustrate, let n = 4. Then the clusters will look like this: 

Vl,V2,V3,V0 Vl,V2,V3,V4 V5 ,3/2,2/3 ,3/4 2/5 ,1/6, J/3, 2/4 2/5,2/6,2/7,2/4 

4 12 3 4 
• • • • • • • • . 



to t\ ti t 3 ti 

In order to include this situation into the setup of Section || (cf . Figure P , we create 
an infinite "caterpillar tree" whose "spine" is formed by the vertices t m , m > 0. 
We thus attach the missing n — 2 "legs" with labels in [n] — {{m — 1) , (m) }, to each 
vertex t m . 

Our next goal is to state conditions on the polynomial F which make it possible 



to assign exchange polynomials satisfying (2.3)-([2.q) to the newly constructed legs. 
The first requirement (cf. ( |2.3| )) is: 

(3.4) The polynomial F is not divisible by any Xj, % £ [n — 1]. 
For m G [n — 1] , we set 

(3.5) Qm = F m |x n <— o = FyXjyi+x , . . . , x n ~i , 0, xi , . . . , x m _i) . 
Our second requirement is 

(3.6) Each Q m is an irreducible element of A[xf 1 , . . . , x^i-J. 

To state our most substantial requirement, we recursively define a sequence of 
polynomials G„_i, . . . , Gx, Go 6 A[xi, . . . , x„_x]; more precisely, each G m will be 
defined up to a multiple in A x . (Later, Gx,... , G„_2 will become the exchange 
polynomials assigned to the "legs" of the caterpillar labeled by n — (0) ; see Fig- 
ure!) 

We set G n _i = F, and obtain each G m _i from G m , as follows. Let 
(3-7) Gm—\ = G m \ Q m . 
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(2) 



<3> 



(0) 



-<o> 



<0) 



<2> 



(0) 



Figure 2. Constructing a caterpillar; n = 4. 



Let L be a Laurent monomial in x\, . . . , x n -\, with coefficient in A, such that 

Cm — 1 



(3.8) 



Gm-l 



is a polynomial in A[xi, . . . , x n -{\ not divisible by any Xi or by any non-invertible 
scalar in A. Such an L is unique up to a multiple in A x . Finally, we set 



(3.9) 



G 



Gm—1 



m—1 



where Q b m is the maximal power of Q m that divides Gm-i- With all this notation, 
our final requirement is: 

(3.10) G = F. 

Theorem 3.1. Let F be a polynomial in the variables xi, . . . ,x n —i with coefficients 



in a unique factorization domain A satisfying conditions ( |3.4| ), (3.6), and (3.10) 
Then every term of the sequence (j/j) defined by the recurrence 

F(y m +i, . . . , y m + n -i) 



Vn 



is a Laurent polynomial in the initial n terms, with coefficients in A. 



Proof. To prove the Laurentness of some y^, we will apply Theorem 2.1 to the 



caterpillar tree constructed as follows. We set thoad = i/v— n+ii this corresponds to 
the first cluster containing y^. As a path from to to thcad, we take a finite segment 



of (|3_3|): 

(3.11) t 



F(i) 



t 2 



<2> 



"(2> 



(N-l) 

F(N-l) 



(N) 



ijV-n+1 ■ 



We then define the exchange polynomial Gj.k-i associated with the leg labeled j 
attached to a vertex tf. on the spine (see Figure ||) by 



Gj,k-i — G{k-j-i)(x{j+i), ■ ■ ■ ,x n ,xi, 



where in the right-hand side, we use the polynomials G\, . . . , G n -i constructed in 

(OhQ above - 



It remains to verify that this exchange pattern satisfies ( |2.3| ), ( p.4[), and (2.5). 
Condition ( [2.3] ) for the edges appearing in (3.11) is immediate from (3.4), while for 

the rest of the edges, it follows from the definition of Gm-i in (3 
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th- 



at) 



FlGURE 3. 



Turning to (2.4), we first note that we may assume i = (0) = n (otherwise apply 
a cyclic shift of indices). Under this assumption, we can identify the polynomials P 
and Qo in (2.4) with the polynomials G m -\ and Q rn i n d3.9| ), for some value of m. 
(The special case of P attached to one of the edges in ( p. 11 ) corresponds to to = 1 , 
and its validity requires (3.10).) Then the condition gcd(G m _i, Q m ) = 1 follows 
from (3.6) a nd the choice of the exponent b in ( 3.9 K_ 



Finally, (2.5) is ensured by the construction (3.7)-(3.9), which was designed 
expressly for t his purpose. A s befo re, the special case of P attached to one of the 
edges in (pll) holds due to ([j~To|) . □ 



In the rest of this section, we give a few applications of Theorem 3.1. In all of 
them, conditions (3.4) and (3.6) are immediate, so we concentrate on the verification 
of (p0|). 



Example 3.2. Let a and b be positive integers, and let the sequence yo,yi,. 
satisfy the recurrence 



!Jk 



Vk-iVk- 



Vk-3 



We claim that every term of the sequence is a Laurent polynomial over Z in yo, yi, 
and y 2 . To prove this, we set n = 3 and construct the polynomials G2, Gi, and Go 
using (FO)-(plj). Initializing G2 = F{x\, x 2 ) — x\x\ + 1, we obtain: 



Q 2 = F(0, Xl ) 

Q 1 = F(x 2 ,Q) 
as desired. 



1, 



Gi=F\ 



Go— Gi 



Q2 — ^1*^2 



1. 



Gi =Gi=x\ + x\, 

G =Go= 1 + x\x\ = F, 



Example 3.3. (Generalized Somos-4 sequence) Let a, b, and c be positive integers, 
and let the sequence yo,yi, . ■ ■ satisfy the recurrence 



Vk 



Vk-M- 



vl-2 



Vk-i 



(The Somos-4 sequence M, introduced by Michael Somos, is the special case a = 
c = 1, b = 2.) Again, each y^ is a Laurent polynomial in the initial terms y , yi, 
j/2, and y 3 . To prove this, we set n = 4 and compute G3, . . . ,Gq using (3.7)— (3.9) 
and beginning with G3 = F = x\x\ + x\: 



Q3 = F(0,xi,x 2 )=x\, 
Q 2 = F(x 3l 0,x 1 ) = xlx1, 
Q 1 = F(x 2 ,x 3 ,0) = x b 3 , 



G, 



x 



2' 



= x a+bc - 
Qa — 3 

r<\ „ _„a+6c , be -b ab+c 



1x2* 



T 2 - 

Gi = 



G2 — X-^~ ~\~ X^X^ 5 
„a_6 1 „ab-\ 



12 3 



Gl I . Qi — 3^1 2 ; 2 a: '3 b "I" ^3 ~*~ > 



G = ^+^=F, 
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and the claim follows. 



Remark 3.4. The Laurent phenomena in Theorems 1.4 and 1.6 can also be proved 



by applying Theorem 3.1: in the former (resp., latter) case, the polynomial F is 
given by F — ax p x n - p + (3x q x n - q + ~/x r x n - r (resp., F = ax p x n -. p + (3x q x n - q ). 
The proofs are straightforward but rather long. Shorter proofs, based on J. Propp's 
idea of viewing one-dimensional recurrences as "projections" of multi-dimensional 
ones, are given in Section ^ below. 

4. TWO- AND THREE-DIMENSIONAL RECURRENCES 

In this section, we use the strategy of Section || to establish the Laurent phenom- 
enon for several recurrences involving two- and three-dimensional arrays. Our first 
example generalizes a construction (and the corresponding Laurentness conjecture) 
suggested by Noam Elkies and communicated by James Propp. Even though the 



Laurent phenomenon in this example can be deduced from Theorem 1.6, we choose 
to give a self-contained treatment, for the sake of exposition. 

Example 4.1. (The knight recurrence) Consider a two-dimensional array {yij)i,j>o 
whose entries satisfy the recurrence 

(4-1) u,,,!/, ; = ayij-iy^j + &yi-\jVi-\j-\ . 

We will prove that every y^ is a Laurent polynomial in the initial entries 

Yhxxt = {yab ■ a < 2 or b < 1}, 
with coefficients in the ring A = Z[a, /?]. 

We will refer to Yi n j t as the initial cluster, even though it is an infinite set. 
Notice, however, that each individual yij only depends on finitely many variables 
{y a b G *init : a<i, b < j}. 

Similarly to Section we will use the exchange relations (Llf) to create a se- 
quence of clusters satisfying the Caterpillar Lemma (Theorem 2.1). 

This is done in the following way. Let us denote by TL = Z> the underlying set 
of indices; for h = G TL, we will write yh = yij ■ The variables of the initial 

cluster have labels in the set 

flinit = {(ij) £ H : «<2or j<l}. 

In Figure [|, the elements of H m n are marked by »'s. 
We introduce the product partial order on TL: 

(4.2) < (i 2 , h) ^> (k < i%) and (ji < j 2 ). 

For an element h = (i,j) e H — H lnlt , let us denote h~ = (i — 2,j — 1); in this 



notation, the exchange relation (4.1 ) expresses the product y^ ■ y^- as a polynomial 
in the variables yh> , for h~ < h' < h. 

We write h~ ~ h, and extend this to an equivalence relation ~ on TL. The 
equivalence class of h is denoted by (h) . These classes are shown as slanted lines in 
Figure ^. All our exchange polynomials will belong to the ring A[a; a :o6H/~]. 

Note that -f/mit has exactly one representative from each equivalence class. We 
will now construct a sequence of subsets Ho = flinit , H\ , H% , . . . , each having this 
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property, using the following recursive rule. Let us fix a particular linear extension 
of the partial order ( |4.2| ) , say, 

d («2, js) (n +ii < «2 + J2) or («! + ji = i 2 +.]2 and i x < i 2 ). 

Restricting this linear ordering to the complement TL — H ln i t of the initial cluster, 
we obtain a numbering of the elements of this complement by positive integers: 

h = (2, 1), hi = (2, 2), h 2 = (3, 1), h 3 = (2, 3), hi = (3, 2), 
h 5 = (4, 1), h 6 = (2, 4), kr = (3, 3), h s = (4, 2), 

and so on. Having constructed H m , we let H m+ \ = H m U {h m } — {/i m }. To 
illustrate, the set Hg is shown in Figure |[ 



1 1 1 i i i i i i i 4 4 4 



FIGURE 5. Indexing set H 9 

We next create the infinite exchange pattern 

(4.3) t — <J^— h — t 2 — t 3 — U • 

f(h ) ^ J (h 1 ) ^(h 2 ) ^(h 3 ) 

(cf. fl3.3|)) The cluster at each point t m is given by x(t m ) = {yh ■ h E H m }; as 
before, each cluster variable yh corresponds to the variable xih)- The exchange 
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polynomial P/m for an edge 



with h = is given by 



(4.4) 



P(h) - ax ((i,j-l)) x ((i-2,j)) +/3x ((l _ lj)) x ((l _ lj _ 1)) . 



Then equations (4.1) become the exchange relations associated with this pattern. 

To establish the Laurent phenomenon, we will complete the caterpillar pattern by 
attaching "legs" to each vertex t m and assigning exch a nge polyn omial s to these legs 
so that the appropriate analogues of conditions ( |3.4[ ), (3.6) and ( |3.10| ) are satisfied. 
Since we now work over the polynomial ring A[x a : a G H/~] in infinitely many 
indeterminates, the number of legs attached to every vertex t m will also be infinite 
(one for every label a different from (/i m _i) and (h m )). This will not matter much 
for our argument though: to prove the Laurentness for any yh m , we will simply 
restrict our attention to the finite part of the infinite caterpillar tree lying between 
to and £hcad = tm+i, and to the legs labeled by (hk) for < k < m. 

The role of conditions (5-4) and (3.6) is now played by the observation that each 
exchange polynomial P/m is not divisible by any variable x a , and furthermore every 
specialization P^)^ is an irreducible element of the Laurent polynomial ring. 



To formulate the analogue of ( 3. 10| ) , let us fix an equivalence class a E H/^ 
and concentrate on defining the exchange polynomials for the legs labeled by a and 
attached to the vertices squeezed between two consecutive occurrences of the label 
a on the spine: 



(4.5) 



Gi a 



G 2 a 



We note that the labels a\, . . . , aAr_i e appearing on the spine between these 

two occurrences of a are distinct. For m = N — 2, N — 3,...,1, we denote by G m 
the exchange polynomial to be associated with the a-labeled leg attached between 
the edges labeled a m and a m +i (cf. (p~5|)). 

The polynomials G m are defined with the help of a recursive procedure analogous 
to (U)-(§T9|^ We initialize G N - 
The step (3.7) is replaced by 



i = P a , and obtain each G m -i from G m , as follows. 



Gm — 1 — G m 



with 
(4.6) 



I cCq.* — ' 



We then compute Gm-i and G r , 



given in the proof of Theorem 3.1 
Laurentness (cf. ( 3 . 1 0| ) ) . 



-i exactly as in 
the equality Go 



(^-(^gl). By the argument 
= P a would imply the desired 
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To simplify computations, we denote the equivalence classes "surrounding" a, as 
shown below: 



(4.7) 



q 


P 


f 


c 


a 


■■■ f 


c 


a 


e 


b ■■■ 


a 


e 


b 


f) 


d ••• 



In other words, if a — then b = — 1)), c = ((i — 1, j)), etc. With this 

notation, we can redraw the pattern (4.5) as follows: 

In \ a q f e d c b a 

(4.8) • • — ^— ■ ■ ■ • — — • • • • • • • • • • • •, 



Gfc-i a 



a 9 Xf + px c x e , 



for appropriate values of k and I. 

We will call a value of m essential if G m -\ ^ G m . We are going to see that the 
essential values of m are those for which a m £ {b, c, e, /}; in the notation of (4.8), 
these values are £ + 1, I, k + 1, and k. 

We initialize Gjv-i = Pa — ax Xf + (3x c x e . The values of m in the interval 
I < m < N are not essential since the variable x am does not enter P a , which is 
furthermore not divisible by Q m (because the latter involves variables absent in P a ). 

The first essential value is m = £ + 1, with a m = b: 

Q i+ i = Pb\ Xa ^o = (ax a Xd + /3x e x g )\ Xa ^o = (3x e x g , 

Ge= P °\ Xb ^^±± - x b 
Gg = aXgXf + Xi,x c . 
Step m = £ (here a m = c) : 

Qi = Pc|x a <-o = (ax e x p + (3x a Xf)\ Xa ^ = ax e x p , 

Gt-1 = Gl\ x ^Q L = OtX g Xf + x h ^f- , 

Gg—i — x c XgXf -\- XfoX e Xp . 

Notice that Gg_\ does not involve Xd, so the value m = k + 2 is not essential, as 
are the rest of the values in the interval k + 1 < m < £. 
Step m = k + 1, with a m = e: 

Qfe+i = f e |x a «-o = (ax c x g + f3x a x b )\ Xa ^ = ax c x g , 

Gk X c XgXj ~t~ -Eb-Ep J 7 

Gfe = x jx e + aXbXp . 
Step m — k, with a m = /: 

Qfc = -P/Ua^o = (ax a x q + /3x c x p )\ Xa ^o = /3x c x p , 

Gfc-1 = /3^c^e + ax b Xf . 
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The values of m in the interval < m < k are not essential since none of the 
corresponding variables x Um appears in Gk-\\ in particular, m — 1 is not essential, 
since Gk-i does not involve x g . Hence 

G = Gk-i = (3x c x e + ax b xf = P a , 

as desired. The Laurentness is proved. 

Remark 4.2. The Laurent phenomenon for the recurrence ( f4.l| ) actually holds in 
greater generality. Specifically, one can replace Ti by any subset of Z 2 which satisfies 



the following analogues of conditions (L2)-(L3) and (1.7)-( |l.8| ): 

(4.9) if h E H, then ti E H whenever h<ti] 

(4.10) for any ti E Ti, the set {h E Ti : h < h'} is finite. 

Then take H init = {h E Ti : h~ £ Ti}. 

The proof of Laurentness only needs one adjustment, concerning the choice of 
a linear extension -< . Specifically, while proving that yt is given by a Laurent 
polynomial, take a finite set Ti^ C Ti containing h and satisfying the conditions 

(4.11) if ti 6 H (h \ then ti' £ H {h) whenever h" < ti and h" e H; 

(4.12) for any ti E H such that ti < h, there exists ti' E Ti {h) such that 
h" > h and h" ~ h. 



(The existence of TL {h) follows from (pT9|) - (f4 . 1 0[) . ) Then define ^ exactly as before 



on the set H {h) ] set ti < h" for any ti E H^and h" eH- H w ; and define d on 
the complement Ti. — Ti^ by an arbitrary linear extension of <. These conditions 
ensure that the sets H m needed in the proof of Laurentness of the given are well 
defined, and that the rest of the proof proceeds smoothly. 

Armed with the techniques developed above in this section, we will now prove 
the main theorems stated in the introduction. 



Proof of Theorem [L2j. Our argument is parallel to that in Example 4.1, so we 
skip the steps which arc identical in both proofs. For simplicity of exposition, we 
present the proof in the special case Ti = Z> ; the case of general Ti requires the 



same adjustments as those described in Remark 4.2 



We define the product partial order < and a compatible linear order ■< on Ti by 

(hJufa) < (*2,J2,&2) ^ (ii < h) and {j x < j 2 ) and (fc x < fc 2 ), 
r< {i<z,h, k 2 ) ^ (ii + ii + h < i 2 + 32 + h) 

or (ii + ji + fci = i 2 + ji + fc 2 and i\ + ji < %2 + .72) 
or {ii + ji = i 2 + j 2 and k\ = k 2 and i\ < i 2 ). 



For h — (i, j, k), we set h = (i — 1, j — l,k — 1); thus, the exchange relation (1.1) 



expresses the product y^ ■ y h - as a polynomial in the variables yf , for h < h' < h. 



All the steps in Example 4.1 leading to the creation of the infinite exchange 



pattern (L3) are repeated verbatim. Instead of (L4), the exchange polynomials 



P(h) along the spine are now given by 

i,k)) 

"{(i-l,j,k)) x ({i,j-l,k-l)) + P x X(i,j-X,h)) x {(i-l,j,k-l)) +l x ((i,j,k-l)) x \(i-l,j-l,k)) ■ 
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The role of (4.7) is now played by Figure which shows the "vicinity" of an 
equivalence class a. This figure displays the orthogonal projection of Ti. along 
the vector (1,1,1). Thus the vertices represent equivalence classes in For 
example, if a — ((i,j,k)}, then 

b={{i,j,k-l)}, c={(i,j-l,k)}, d=((i-l,j,k)), 

e={(i,j-l,k-l)), f= {(i-l,j,k-l)}, g=((i-lj-l,k)}. 

With this notation, we have: 

P a = aXdX e + f3x c Xf + JX b X g . 




10 

Figure 6. The cube recurrence 



With the polynomials Gi,G2, ■ ■ ■ defined as in (4.5), the essential values of m 
are now those for which a m G {&, c, d, e, /, <?}. (The verification that the rest of the 
values are not essential is left to the reader.) We denote these values by mi, . . . , mg , 
respectively. 

The computation of the polynomials G m begins by initializing 

Gjv-i = P a = ax d x e + (3x c x f + jx b x g . 
Step m = mi, a m = b: 

Qva\ — Pb\x a < — — OLX^Xq -\- f3x e Xp , 

G mi -i= G mi\ Xbt Qm! = ax d x e + [3x c x f + 7 — -x g ; 

G mi _i = ax b x d x e + (3x b x c Xf + ajx f x g x q + (3jx e x g x p . 
Step m = TO2, Cbm = c: 

Qm2 Pc\x a < — OLXgXf -\- 'yXeXs, 

Gm 2 -i= ax b x d x e + f3x b aXgXr ^ x ' Xs Xf + ajx f x g x q + (3"/x e x g x p ; 

+(3^x b x e x fX s + ajx c x fx g x q +f3^fx C X e XgXp . 
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Step m = 7713, a m = d: 

/3X g Xy+JXfX u 

(-Tm.3 — 1 — QLX\>X C ^ X e 

+a/3xbX fXgX r + PjXbX e XfX s + ajx c x fXgX q + /3jx c x e x g Xp ; 

x e XgX v -\- oi'yxi)X c x € _XjXii H- f3 r yxbXdX e x j x s ~\~ i3 r )x c X(ix e .XgXp 
+af3xbXdXfX g x r + a~/x c XdX fXgX q . 
Step m = m4, a m = e: 

Qm^ ' -^ela^a^ — ~~ $X\yX r -\- ^fX c Xq , 

G m4 -i= -^{a(3x b x c x g x v + a-yx b x c XfX u + (3-fx b x d XfX s + f3jx c x d x g x p ) 

-\~OiXdX fXgQrri4 i 
G 1114—1 — ajXbXcXfXu + /3jXbXdXfX 8 + aXdXeXfXg 
I Oi. ^3 X I) X (• X g X d I X qX dX gX rp • 

Step m = m 5l a m = /: 

Q ms = Pf\ Xa ^ = axbx v + jXdXp ; 

Cms— 1 — x f [^^f'^b'^C'Eu H~ P"^fXf}XdX s ~\~ QtXdX e Xg^ H~ ^c-^gQms 7 

G m5 _l = ax d x e Xg + f3x c x f Xg + cryxfcXcO^. + (3^x b x d x s . 
Step m = me, a m = 5: 

Qmg — ■^fl'l^a* - — Q;X C X W H - {3XdX s , 
Gm 6 -1= —^(ctXdXe + (3x c Xf) + JXbQmz \ 



G m6 -i = ax d x e + [3x c x f + jx b x g = P a , 
completing the proof. □ 



We will now deduce the Gale- Robinson conjecture from Theorem 1.2. 

Proof of Theorem |1.4| . To prove the Laurentness of a given element of 
the Gale-Robinson sequence (j/ m ), we define the array {zijk)(i : j.k)£H by setting 
Zijk = VN+pi+qj+rk , with the indexing set 

H = H{N) = {(i,j,k) G Z 3 : N + pi + qj + rk > 0} . 



Then (1.4) implies that the z^ satisfy the cube recurrence (1.1). Note that H 
satisfies the conditions ( |1.2| )-( fL3| ). Thus Theorem 1.2 applies to (zyfc), with H init — 
{(a, b, c) G Z 3 : < N + pa + qb + rc < n}. It remains to note that yjv = ^ooo , 
while for any (a, 6, c) G -ffinit , we have z a b c = y m with < m < n. □ 



Proof of Theorem 1.6. This theorem is proved by the same argument as The- 



orem |L2] . We treat the Mills-Robbins-Rumsey special case fll.9Q (cf. also (1.6)); 
similarly to Theorem |1.2| , the case of general 7i requires the standard adjustments 
described in Remark [4. 2| . We use the partial order on the lattice L defined by 

(i,j,k) < (i',j',k') : \i' + \f -j\ < k' ~ k . 
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For h — k) E L, we set h~ = (i, j, k — 2), and define the equivalence relation ~ 
accordingly. Figure [?] shows equivalence classes "surrounding" a given class a (cf. 
Figure ||) . 



Phi 



I i 

d 


! c 

a 


(! i 

S 

I ( 


> ( 

e 

9 c 



Figure 7. 



The initialization polynomial Gn-i — P a is given by P a = ax c Xd + f3xi,x e . The 
table below displays a m , Q m , Gm-i, and G m _i for all essential values of m. 

Qrn Qm Cm — 1 G m _i 



C 



QtXpX q (XX qX d ~\~ Oift Xg XfjX qX d + Px e X p Xg 

j3XqX r ft — X{)X(i -\- f3x e XpXq XfoX dX r -\- X c X e Xp 

ftXpX s ft ~ XfjX r H - X c X e Xp ftXfjX r X s ~\~ X c Xd%e 

^ (XtXjy X" g ft X/ iX X^ q I ~ Xj q Xj d ft \) *^ £ I CXtjC Xt d 

We see that Go = G e _i = P a , completing the proof. □ 
Proof of Theorem |1.8|. The proof mimics the above proof of Theorem 11.41. To 



prove the Laurentness of an element i/n of the sequence (y m ) satisfying (1.10), 
we define the array (Zijk)(i,j,k)eH by setting z i]k = y N +e(ij,k) , where £(i,j,k) = 



pi — qj . The indexing set H is now given by 



H = H(N) = {{i,j,k)eZ 3 : N + t{i, j, k) > 0} . 

Then ( l.lOj ) implies that the Zijk satisfy t he o ct ahe dron recurrence (1.5). It is easy 
to check that Tt satisfies the conditions (1.7)- (|l.8[) . Thus Theorem F6 applies to 
(zijk), with ifinit = {(a, b, c) E L : < N + £(a, b, c) < n}, and the theorem follows. 
□ 

We conclude this section by a couple of examples in which the Laurent phenom- 
enon is established by the same technique as above. In each case, we provide: 

• a picture of the equivalence classes "surrounding" a given class a, which plays 



the role of (4.7) in Example 4T 



• the initialization polynomial Gtv_i = P a : 

• a table showing a m , Q m , Gm-i, and G m -i for all essential values of m. 

Example 4.3. (Frieze patterns) The generalized frieze pattern recurrence (cf., e.g., 
§ 0) is 
(4.13) 
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where s £ {1,-1}. To prove Laurentness (over Z [/?]), refer to Figure ||. Then 
P a = e xjy x c + [3, and the essential steps are: 

dm Qm Gm—1 G m —\ 



b p £|£k +/ 3 ex c + x b 
c (1 Z£+Xb (i + e~ x x h x c 




Figure 8. 



Example 4.4. (Number walls) Consider the 2-dimensional recurrence 

( 4 - 14 ) yijVij-2 = yf^^yl+Lj-! + , 

where p, q, and r are nonnegative integers. To prove Laurentness, refer to Figure 
Then P a = x v d x\ + x q c , and the essential steps are: 

Q>m Qm Gm — 1 G m — i 



r Q 



p qr 
X d X f 




Figure 9. 



Remark 4.5. As pointed out by J. Propp, the Laurent phenomenon for certain 
special cases of Examples 4.3 and 4.4 can be obtained by specialization of Exam- 
ple pi 
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5. Homogeneous exchange patterns 



In this section, we deduce Theorem 1.10 and a number of similar results from 



the following corollary of Theorem 2.1 



Corollary 5.1. Let A be a unique factorization domain. Assume that a collection 
of nonzero polynomials P\, . . . , P n G A[xi, . . . , x n ] satisfies the following conditions: 

(5.1) Each Pk does not depend on x k , and is not divisible by any x\, i £ [n]. 

(5.2) For any i ^ j, the polynomials Pji d = (Pj)\ Xi =o and Pi are coprime. 

(5.3) For any i ^ j, we have 

L ■ Pji ■ Pi — Pi\ p }i , 

where b is a nonnegative integer, and L is a Laurent monomial whose 
coefficient lies in A and is coprime with Pi . 
Let us define the rational transformations Fi, i € [n], by 

Pi 

Pi ■ (^1) • * • j ) 1 * (^l) • ■ ■ ) 1: ) ^i+lj ■ ■ • ! ^n) ■ 

Then any composition of the form Fi 1 o • • ■ o Fi m is given by Laurent polynomials 
with coefficients in A. 

Proof. Let T„ denote a regular tree of degree n whose edges are labeled by 
elements of [n] so that all edges incident to a given vertex have different labels. 
Assigning Pi as an exchange polynomial for every edge of T„ labeled by i, we 



obtain a "homogeneous" exchange pattern on T„ satisfying conditions ( |2.3| )-(2.5) 
in Theorem 2.1. This implies the desired Laurentness. □ 

Example 5.2. Let n > 3 be an integer, and let P be a quadratic form given by 

P(x\, . . . , Xyi} X^ ~\~ ' • • ~\~ X n ~~\~ ^ (XijXiX j . 

i<j 

Theorem [LTo| is a special case of Corollary [Tl] for Pj = P| and A = Z[a.y : i < 

j]. Conditions (pTl|)-(5.2) are clear. To verify ([5~3]), note that 

Pi = Pji + X 2 j + X.j a kjXk + a 3l X l) ' 

k t 

where k (resp. £) runs over all indices such that k ^ i and k < j (resp. £ ^ i and 
£ > j). It follows that 



p *\ x ,_2i= p j* + -3 + ^r(Yl a ^k + a H x l) = ^ p 



X j X 3 K k e X j 



verifying (|5.3|) 



In the remainder of this section, we list a few more applications of Corollary 5.1 
In each case, the verification of its conditions is straightforward. 

Example 5.3. Let P and Q be monic palindromic polynomials in one variable: 
P(x) = (1 + x d ) + ai(x + x d - r ) + a 2 {x 2 + x d - 2 ) + ... : 

Q{x) = (1 + x e ) + [3i(x + X s - 1 ) + (3 2 {x 2 + x e - 2 ) + ... . 
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Then every member of the sequence yo,yi, ■ ■ ■ defined by the recurrence 

> 2 P(y fc -!/A) .... 
it k is odd; 

= , Vk-2 

< X iQjvH-i/ri .... 

it fc is even 

Vk-2 

is a Laurent polynomial in yo and z/i with coefficients in A — Z[A ±1 ,/i ±1 , aj,/3j]. 
This follows from Corollary 5.1 with n = 2, Pi = /x 2 P(x 2 /A), and P2 = X 2 Q(xi/ fi). 

Example 5.4. Consider the sequence yo,y±,... defined by the recurrence 

y 2 k -! + cy k -i + d 

(5.4) y k = . 

Vk-2 

Every term of this sequence is a Laurent polynomial in y and yi with coefficients 
in Z[c, d]. 

Example 5.5. Define the rational transformations Pi,P 2 ,F3 by 
rp t \ , x 2 + x^ + x 2 x 3 

i<i : (xi,x 2 ,x 3 ) ( , x 2 , x 3 ), 

Xl 

(5.5) F 2 : fax, X2,x 3 ) h-> ( xi, + Xa , x 3 ), 
tp 1 \ ( x 2 + x\ + xlxi 

F 3 ■ (Xl,X 2 ,X 3 ) ( Xl, x 2 , ). 

Then any composition P ix o P i2 o • • • is given by (xi, x 2 , X3) 1— > (Gi, G 2 , G3), where 
Gi, G 2 ,G3 are Laurent polynomials in xi,x 2 ,X3 over Z. 
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